Topography and instability of monolayers near domain boundaries 



O 

o 

(N 



O 



> 

m 

00 



I 



X 



H. Diamant/ T. A. Witten,i C. Ege,^ A. Gopal,^ and K. Y. C. Lee^ 
^ James Franck Institute and Department of Physics, University of Chicago, 
5640 South Ellis Avenue, Chicago, IL 60637 
^Department of Chemistry and Institute for Biophysical Dynamics, University of Chicago, 
5735 South Ellis Avenue, Chicago, IL 60637 
(February 9, 2001) 



We theoretically study the topography of a biphasic surfactant monolayer in the vicinity of domain 
boundaries. The differing elastic properties of the two phases generally lead to a nonflat topography 
of "mesas" , where domains of one phase are elevated with respect to the other phase. The mesas are 
' steep but low, having heights of up to 10 nm. As the monolayer is laterally compressed, the mesas 

-jr , develop overhangs and eventually become unstable at a surface tension of about K{5cQf ('^cq being 

the difference in spontaneous curvature and K a bending modulus). In addition, the boundary is 
^ I ^ ' found to undergo a topography-induced rippling instability upon compression, if its line tension is 

smaller than about KScq. The effect of diffuse boundaries on these features and the topographic 
behavior near a critical point are also examined. We discuss the relevance of our findings to several 
experimental observations related to surfactant monolayers: (i) small topographic features recently 
found near domain boundaries; (ii) folding behavior observed in mixed phospholipid monolayers and 
55 ' model lung surfactants; (iii) roughening of domain boundaries seen under lateral compression; (iv) 

the absence of biphasic structures in tensionless surfactant films. 

68.18.-g,82.70.Uv,87.68.+z 
O ' I. INTRODUCTION 
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Monolayers of amphiphilic molecules (surfactants) at water/air or water/oil interfaces are used in numerous appli- 
cations to reduce interfacial tension, control wetting properties, stabilize emulsions and foams, etc. Monolayers 
of biological surfactants (phospholipids) are commonly studied as models for the surfaces of cell membranes and are 
also encountered in various biological systems An important example is the lung surfactant monolayer covering 
the alveoli in lungs, whose main function is to lower the surface tension of the lungs, thereby drastically reducing the 
mechanical work required for breathing |4| . 

Amphiphilic monolayers generally have a finite spontaneous curvature arising from the asymmetry of the molecules 
as well as the asymmetry with respect to electrostatic interactions {i.e., the differing dielectric properties of the 
polar and nonpolar phases forming the interface) Q. Despite this tendency to bend, homogeneous monolayers are 
almost always flattened by the water/air or water/oil interfacial tension. Only for very low (sometimes even negative 
1^) tension does a nonflat conformation become energetically favorable for a homogeneous monolayer. (This is 
achieved, e.g., by extensive lateral compression.) Such a reversible departure from a flat, two-dimensional state to a 
. three-dimensional conformation is referred to as the buckling transition and has drawn considerable attention [^|-^. 
(~| ' However, it is not commonly observed in practice |7j, since it is usually preceded by other modes of collapse such as 
Q , monolayer breakage into multilayers HJl^ and ejection of vesicles or micelles The possibility to explore the third 
O ' dimension upon compression is of particular interest in the case of lung surfactant monolayers, which are required to 
^ , change their projected area significantly during the compression-expansion cycle of breathing. 

The two-dimensional fluid comprising a monolayer may separate into domains of different coexisting phases. Single- 
component monolayers exhibit gas/liquid-expanded, liquid-expanded/liquid-condensed and liquid-condensed/solid co- 
existence whereas mixed monolayers may form domains of differing composition. A special, well-studied property 
of surfactant monolayers is the stabilization of finite domains and modulated phases due to long-range electrostatic 
interactions p2|Jl^ ]. The coupling between lateral variations in composition and curvature was thoroughly studied 
as well [p|,p^[-pO|7rna-iiily with regard to various domain structures on surfaces and shape transformations of bilayer 
vesicles. 

Despite extensive research on surfactant monolayers there are important features, in particular of biphasic mono- 
layers, which are not well understood. Recent experiments on mixed phospholipid monolayers have revealed a new 



type of local folding upon compression, which is believed to be important for the function of lungs |21| (see Fig. |^). 
Another observed feature is the appearance of rough domain boundaries upon compression p2| (see Fig. P). 

In the current work we study the relation between lateral domain structure and monolayer topography in more 
detail, focusing on the conformation of monolayers in the vicinity of domain boundaries |23||. We thereby try to 
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shed some light on the unexplained features mentioned above. Domains of different density and/or composition in a 
biphasic monolayer should generally have differing elastic properties. The requirements of mechanical equilibrium and 
smoothness of the monolayer surface lead to nonflat conformations attempting to "reconcile" the different properties 
of the contiguous domains. These simple mechanical considerations result in a surprisingly rich behavior, including 
the formation of overhangs and emergence of instabilities, as discussed in the following sections. 

The Monge representation and linearization of profile equations have been ubiquitously used in the theoretical 
modeling of monolayers and membranes [^j . These mathematical simplifications describe the topography by a single- 
valued height function assumed to have moderate slopes. By contrast, the phenomena discussed in the current work 
involve, in an essential way, steep slopes and overhangs. We thus avoid the Monge representation and solve the 
nonlinear profile equations. In order for the mathematics to remain tractable we resort to another simplification — the 
profile is assumed to be uniform along one lateral direction (namely, the direction parallel to the domain boundary) — 
rendering the nonlinear equations one-dimensional. This constraint is further discussed in the next section; it is 
partially relaxed in the treatment of boundary rippling in Sec.^^ 

The basic model and its assumptions are presented in Sec. We then proceed in Sec. Ill to review the simplest 
case of an infinitely sharp, straight domain boundary |23| . Despite its simplicity, this limiting case demonstrates most 
of our qualitative results. The calculation is refined in Sec. where a boundary of nonzero thickness is considered. 
In Sec. 1^ we study the stability of a straight domain boundary to lateral und ulations. The topographic behavior of a 



monolayer near its critical point is examined in Sec. |VI|. Finally, in Sec. |VI|, we discuss the various results and their 
relevance to experiments. 



II. THE MODEL 



Four length scales are distinguished when studying the elasticity of a biphasic monolayer: the typical domain size, 
L, the width of a domain boundary, 5, the typical spontaneous radius of curvature, c^^, and the elastic length, 
A — (Kj^Yl'^ , determining the lateral extent of height variations [K being the bending modulus and 7 the surface 
tension). An important observation is that in most practical cases the domain size is much larger than all other length 
scales — L is typically of order 10 /xm, whereas ^, c^^ and A are of order 1-10 nm. This allows us to focus on a single, 
straight boundary between two large domains and regard the centers of the domains as infinitely far away. We thus 
represent the boundary region as a surface whose far left and far right have different spontaneous curvatures, cqi and 
C02 I p^ . (Throughout this paper we assume, without loss of generality, Cqi > C02.) The surface is uniform in the y 
direction parallel to the boundary but can curve in the perpendicular x direction, as depicted in Fig. ^A. 

In fact, an inflected conformation perpendicular to the boundary is a straightforward consequence of such a lateral 
structure. Far away from the boundary (i.e., at the centers of the two contiguous domains) the surface is flat. Because 
of the nonzero spontaneous curvature the far left side is subjected to a bending moment of Kcq\, supplemented by a 
tensile torque of — 7/11, where — h\ is the height of this side relative to the boundary (see Fig.0B). Similarly, a moment 
of Kcqi ~\- 7/12 is exerted on the far right side. Mechanical equilibrium requires that these two moments balance each 
other, i.e., 

h = hi+h2^ KSco/j = X^Sco, (1) 

where Scq = cpi — Cq2. Thus, an inflected conformation with a finite height difference occurs for any finite Scq and 7. 
A measure of the inflection sharpness is h/X = XScq. As compression increases {i.e., 7 decreases ^^), the inflection 
becomes higher and sharper. Note that the integrated height difference is insensitive to details of the inner boundary 
region. Equatio n ([l| ) can therefore serve as a rigorous sum rule for more detailed calculations such as those presented 
in Sees, g and ^ 

The elastic energy to be minimized in order to find the monolayer topography pq] is 



G^J dA Qi<rc2 - cqc^ + 7 y d{A-Ap) + Tj d{R 



V). (2) 



where A denotes the monolayer surface, A-p its projection onto the xy reference plane, c the local surface curvature, 
R the trajectory of the domain boundary, and r the line tension of the boundary. The functional (^ has been defined 
such that a flat surface with a straight domain boundary running parallel to the y axis is the reference state of zero 



energy. Since the surface we consider in Sees. Ill and IV is uniform in the y direction, we may represent it by the 
local angle 9{s) it makes with the reference xy plane at curvilinear distance s from the boundary (see Fig. |i|b). The 
energy functional (0) is then rewritten as 
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G 
L 



K0^ - KcqO + 7(1 — cos ( 



(3) 



where L is the length of the boundary, a dot denotes d/ds, and cq is not regarded as a constant but varies with s. 
In order for the surface to be smooth everywhere we require that 9{s) be a continuous function [^7| . It is useful to 
notice that, despite the linear term in 6, the functional (H) is invariant under s-reversal, s — s. This is because 
we can break any plausible choice of co(s) into a constant, (cqi + co2)/2, plus an odd function of s and, assuming 
that the surface becomes flat far away from the boundary, the constant term does not contribute to the integral. 
In other words, we may specialize to the case Cqi = — co2- Hence, minimization will necessarily produce symmetric 
angle profiles 9{s) and antisymmetric topographies h{s). (Note that this argument does not hold if K is taken as 
nonuniform p3|.) 

III. SHARP DOMAIN BOUNDARY 

We begin with the simple case where the boundary thickness ^ is much smaller than the other length scales, 
extending our earlier work psf . (We show in Sec. ^ that the condition for this sharp limit is in fact weaker, the 
requirement being merely ^Scq <^ 1.) In this limit the boundary may be regarded as infinitely sharp, accompanied by 
a step-function jump in spontaneous curvature. 



Co 



Cqi s < 

C02 s > 0. 



Substituting this spatial variation of Cq in Eq. 

as 



-ms)] 



ds 



we can integrate the linear term in 9. The energy is then rewritten 
1 



'-K9^ +"f{l-cose) 



K6cn9, 



(4) 



where 9q = 9(s = 0) is the maximum inflection angle. 

The integral in Eq. (Q) has the familiar form of the Sine-Gordon action. Variation with respect to 9{s ^ 0) gives a 
Sine-Gordon profile equation. 



First integration of Eq. (|^) yields 



= A"^sin( 



2X-^sm{9/2) s<0 
-2X-^sm{9/2) s>0. 



(5) 



(6) 



In the current, simple case, due to the boundary conditions 9{s —>■ ±00) = 0, second integration can be carried out as 
well. The following soliton profile is obtained: 



tan - = tan — exp(— |s|/A). 



(7) 



Finally, we need a condition for the jump in curvature at the sharp boundary. This is derived mathematically by 
either integrating 9 of Eq. (|^) over an infinitesimal length around the boundary, or taking the variation of g with 
respect to 9o- The same result is obtained, nonetheless, by a simple moment-balance argument: the bending moment 
acting on the boundary from the left, K[9{0^) — cqi], must balance the one acting from the right, K[9{0^) — C02]. 
Hence, 



9(0-) - 9(0+) = 6co, 



which determines 9c) as 



X5co 



h 



(8) 



The dependence of the slope on the ratio h/X in the sharp-boundary limit readily follows from dimensional analysis. 
Yet, the exact nonlinear dependence given by Eq. (||) is essential for our main results, as will be demonstrated below. 
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Thus, we infer again that for any finite Scq and 7 the monolayer attains an inflected shape whose profile is given 
by Eqs. (0) and (^. Integrating J^^dssm6{s), one verifies that the general sum rule for the total height difference, 
Eq. (|l|), is satisfied. The energy of the inflected conformation is calculated by substituting the obtained proflle back 
in Eq. (§, 

g{9o)^KSco(2t'^n^-eo], (9) 



which, as expected, is negative for 60 < tt, i.e., the inflected shape is favored over the flat one. We also calculate the 
projected area as a function of tension, either by integration, Ap = L dscos0, or by the following derivative, 



and the lateral compressibility. 



AT ^-^P ^9 32 . 00 . 2 ^0 



dAL ^ 1 [4sin(V2)]3sin(3go/2) 
97 K{6co)^ cos(e'o/2) ■ ^ ' 



As the monolayer is compressed by progressively decreasing 7 or increasing AL (depending on the experimental 
setup), the inflected profile becomes sharper (larger 6*0, Eq. (H)), higher (larger h, Eq. (|l|)) and more favorable (lower 
g, Eq. (||)). The process is demonstrated in Fig. g. For 7 < ^K{6co)^ the inflection angle becomes larger than tt/2 
and a stable overhang forms. However, beyond a critical value of compression, 

7 < 7c = AL>ALe-^, (12) 

Id oco 

Eq. (^ has no solution. This occurs in the current calculation when ~ ^Oc = tt, at which point the inflection 
remains finite, h — he — 16/i5co, yet the lateral compressibility diverges, 

C^(7r-0o)-'^(7-7c)-^/^ (13) 

Note that this instability is revealed only in the nonlinear representation of the elastic problem {e.g., via Eq. (||) which 
is nonlinear in ^o)- A theory relying on the Monge representation and moderate slopes would inevitably miss it. 

The divergence of lateral compressibility implies that extra surface area can be pulled into the boundary region 
without resistance. Thus, one expects the monolayer to attain a folded structure that will be stabilized by higher, 
nonlinear-elastic terms. A detailed description of this folded shape is beyond the scope of the current work. Moreover, 
the transition from an inflected to a folded state is treated here as a spinodal-like instability. A more detailed 
study might yield a folded structure as energetically favorable prior to this instability, i.e., a binodal-like transition 
preceding the one treated here. In a macroscopic measurement of a pressure-area isotherm the instability will appear 
as a plateau, much like any other Ist-order transition. Hence, distinguishing it from other mechanisms requires 
microscopic visualization of the monolayer [E8|. 



IV. DIFFUSE DOMAIN BOUNDARY 



We now examine the effect of a domain boundary of finite thickness. That is, we suppose that rather than changing 
abruptly, cg varies gradually from cqi to C02 over a distance ^. We begin with a heuristic argument for the limit of 
a very diffuse boundary, where ^ is much larger than A and {Sco)~^. Consider a small element ds in the boundary 
region. The difference in bending moments acting on its two ends, Kd{9 — cq) = K{6 — co)ds, is balanced by a 
difference in tensile torque jdh = jsmOds. In the limit of very large ^ the curvature changes very slowly, such that 
6 ^ Oq/S,"^ is negligible compared to cq ~ —5co/£,. We thus obtain 

sin6lo - A^^co/C = V6 (14) 

The mesa slope in this diffuse limit, unsurprisingly, depends on ft,/^ rather than h/\. We can infer from Eq. ( p^ ) 
several less expected results as well. The equation has no solution for > ^Oc — t^/2 (rather than tt), whereupon 
the mesa height is of the order of ^, i.e., very large compared to A. Thus, at the critical compression the entire wide 
boundary lies almost vertically, having a small lateral dimension. The critical tension is 
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7c ~ K5cq/^, 



(15) 



i.e., much smaller than its value in the sharp limit (~ Kidc^Y'), implying that the monolayer can sustain much 
higher compression than in the sharp case. (Naturally, as ^ tends to infinity one expects the resulting almost-uniform 
monolayer to become increasingly stable to the heterogeneity-driven folding.) 

In order to check these results in more detail and study the crossover between the sharp and diffuse limits, we now 
turn to a detailed treatment of a simple, specific example. Let us assume that the spontaneous curvature cq changes 
linearly across the boundary. 



Co 



Substituting the spatial dependence of cq in the energy functional 
we obtain the following profile equations: 



s < -e/2 
fco(s/e + l/2) -e/2<s<e/2 

s > e/2. 

and taking the variation with respect to 0{s), 



X ^ sin ( 
A^^ sin( 



Sco/^ 



M > e/2 

isi <e/2. 



(16) 



The boundary conditions in the current case are continuity of 9 and 9 at s = ±^/2 and, as before, flatness at infinity. 
First integration of Eq. (|l^) gives 



n2 



AX- 
AX 



'2 ■ 2 

^ sm 
-2sin2 



(2<5co/O(0- 



|s|>^/2 
h) <?/2, 



(17) 



where 9\ = 9{—(^/2) = 9(^/2). As in Sec. Ill, we can now calculate the energy of the inflected conformation as a 
function of 9o and 9i by substituting the profile (^) back in the energy functional (H). The result is 



KSco 



XScq 









4 sin^ — + 


/ d9 


sin^ 


I 4 







2e 



1/2^ 



(18) 



Second integration, which would yield the topographic profile, cannot be analytically performed in the current case. 
Instead, we seek the equation for 9o which replaces Eq. (||) for finding the point of instability. Minimizing g of Eq. ( p^ ) 
with respect to 9q gives one relation between 9q and 9i, 



2? . 2 ^0 



This equation reflects a moment balance for the section — oo < s < 0; it is also obtained by setting ^(0) 
Eq. ( p^ ) . Minimization with respect to 9i yields the second relation between Oq and 9i , 



00 r 



2^ .2^ 
A^'^" 2 



-1/2 



d9 



(19) 
in 

(20) 



It stems from the moment balance on the section — ^/2 < s < and can be also obtained from integration of the 
profile equation ( |l6| ) along this section. 

Equations (p9|) and (po[) can be solved numerically for ^o- As in Sec. Ill, one finds a minimum value j = -f^ {i.e., 



a maximum value of A)beyond which there is no solution to the equations and the monolayer becomes unstable. 
From the fact that the various parameters are grouped in Eqs. (^|) and ( pO| ) into two dimensionless terms, ^Scq and 
^/(A^(5co) = £,/h, it follows that the critical tension must satisfy the scaling law 



7c 



KSco 



Fiasco) 



(21) 



where F{x) is a certain "universal" , dimensionless function. In Fig. ||A we have plotted the function F{x) as obtained 
from numerical integration of Eqs. ( p^ ) and (pO[). 

In the sharp limit, S,Scq 1, 7c must become independent of ^. Hence, F(x) is linear for small x so as to get 
7c ^ K{6co)^. Indeed, in this limit Eq. (EG) reduces to 9o — 9i ~ which, together with Eq. (Il9), recovers 
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the results of Sec. |l| (cf. Eqs. (|) and sin(6lo/2) ~ XSco/A = h/{4X), Oqc ^ tt, and 7c ~ -^K{ScQf. We thus 

conclude that the results of Sec. HI are vahd as long as ^ ^ (Scq)^^ . 



In the diffuse limit, ^(5co » 1, Eq. ( |2C| ) reduces to 

sin6io ~ A^fco/e = VC, (22) 

as was anticipated in Eq. This leads to Oqc ~ 7r/2, jc — K5cq/^. Hence F{x) ~ 1 for x 3> 1, as is verified in 

Fig. |A. 

To summarize these results: 

K5cq 
7c = F(f(5co) 

Recall that we have obtained the results in both the sharp and diffuse limits independently of the detailed shape of 
co(s). Hence, varying the spatial dependence of cq would merely affect the exact shape of F{x) in between these 
limits. The more diffuse the boundary, the higher the compression required for folding, the higher the mesa wall at 
the instability, and the smaller the critical inflection angle. Diffuse boundaries thus allow a biphasic monolayer to 
withstand stronger compression and higher mesas. Despite the smaller inflection angle there is always an overhang 
topography at the instability, i.e., 7r/2 < ^oc < i"- The dependence of the critical inflection angle on boundary 
thickness is shown in Fig. pp3. 



V. INSTABILITY OF THE DOMAIN BOUNDARY 



Until now we have considered only topographies which do not vary along the direction of the domain boundary, 
and thus do not affect its length. The departure from a flat conformation near a domain boundary, as studied in the 
previous sections, is energetically beneficial, i.e., the inflection energy per unit length, g, is negative. Hence, as far 
as the topographic effect is concerned, it would be favorable to increase the boundary length. In other words, the 
topography effectively reduces the line tension of the phase boundary, the reduction being given by g of Eqs. or 
( p^ for a sharp or diffuse boundary, respectively. Consequently, if the bare line tension of the boundary |^, r, is 
smaller than the maximum value of \g\, then, for a certain inflection angle 9q < Oqc, the effective line tension will turn 
negative and one expects the boundary to ripple. Assuming hereafter a sharp boundary, we obtain the condition for 
rippling by setting 6*0 = 6*00 = tt in Eq. (|) Q, 

T < Tc = (tt - 2)KdcQ. (24) 

The inflection angle and surface tension required for rippling, ^or and 7r, are then obtained from the equations 



70r 



2 tan — 



4 KSco 



sm {00,/ 2) 

where, as defined in Eq. ([l^), 7c = j^K{Scq)'^ . The diagram in Fig. |^ summarizes the results concerning the 
topographic transitions near a sharp boundary as a function of surface tension and line tension. 



Let us now examine the spatial form of the rippling transition. In Sec. HI we assumed a straight, sharp domain 
boundary, which can be represented in Cartesian coordinates as the line R(y) = {x = 0,y, z = 0). We now wish 
to perturb the inflected conformation by considering a boundary that slightly wiggles with amplitude a and wave 
number q. The full three-dimensional problem is formidable. We therefore restrict ourselves to a simple subset of 
perturbations — uniform displacements of the inflected shape in the x direction, whose magnitude undulates in the 
y direction (see Fig. Since we do not exhaust all available conformations, the minimum energy that we are 
about to calculate might be higher than the true minimum. Hence, the following results should be considered as an 
upper-bound estimate for the rippling instability. Nevertheless, this estimate is expected to be good as long as the 
wiggling wavelength is much larger than the inflection extent, qX <^ 1. In this limit the two lateral length scales can 



be separated, as has been done in Sec. IH, and one expects the preferred perturbations to resemble that of Fig. 
Employing this simplification, we can represent the perturbed boundary by the curve R(y) = {x = a sin qy,y, z — 0), 
and conveniently parameterize the monolayer surface as 
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r(s,t) = {x{s,t),y{s,t),z{s,t)) 

x{s,t) = / cos 6 {s' ) ds' + a sin qt, y{s,t) = t, z{s,t) 



sin9{s')ds', 



such that in the "material coordinates" (s, t) the boundary hne is again given by R(t) = {s — 0,t, z = 0). 

In order to use the energy functional (0) we need to represent the geometrical parameters of the surface — its mean 
curvature c(s,t), area element dA{s,t), projected area element dAp(s,t), and boundary arc-length dR{s,t) — using 
the new coordinates. This technical calculation is presented in the Appendix. 

Substituting Eqs. (A3)-(A6) into the elastic energy expression, Eq. (||), and expanding to second order in the 
rippling amplitude a, we obtain 



9[9{s)] EE G/L = 5^°) + a2g(i) 



(26) 



where g^'^^9{s)] is the energy functional in the straight case, given by Eq. (H), and L = J dt. The coefficients r' and 
b are functional of the topography 9{s): 



K 2~ - sin^ 



f - 2Kco cos^ I 



7 sm 



b = K dssin^6l. 



(27) 

They act as effective line tension and bending modulus, respectively. The one-dimensional bending modulus b is 
proportional to K, the two-dimensional modulus of the sheet; when the boundary curves the monolayer must bend 
with it (see Fig. leading to a cost in bending energy. 

In principle one could now minimize g of Eq. ( P6[ ) with respect to 9{s) and find the perturbed shape, 9{s) = 
6i("'(s) + a26l(i)(s)- However, since Sg^^^ /S9'-^^ — 0, substituting the perturbed shape back in g would yield, up to 
order , g = g^*^") [9^^^] -l-a^g^^^ [0*^°^]. Thus, if we are merely interested in the perturbed energy, we may just substitute 
in Eq. ( p7[ ) the unperturbed topography 0(°)(s) as found in Sec. II] (Eqs. (H) and (^)). This yields 



r' = T - K6co 9o - 2 tan 



, 32K . 6^0 



cos 



70 

4 

3 ^0 



(28) 



When the effective line tension vanishes, r' = 0, there is a 5 = {i.e., 2nd-order) rippling transition, as already 
anticipated in Eqs. (|24) and (Eq). (Strictly speaking, since the domain boundary is finite and closed, the transition 



is encountered only at t' — ^7?6/i^, i.e., for the lowest-order undulation of g = '^/L.) The rippling of the one- 
dimensional boundary is thus analogous to the Euler buckling of an elastic rod | |3l[ | . 

Upon further compression, or if the monolayer is "quenched" to r' < 0, all modes satisfying 



q<q 



-T 



1/2 



aV2 



Sco 



9o - 2tan(6>o/4) - t/{KSco) 
sin(0o/2)[l-cos3(^?o/2)] 



1/2 



(29) 



become unstable and their amplitudes start growing. We expect the observed unstable modes to have roughly the 
same scale as the upper bound q*. The scale of the rippling wavelength is thus set by {Sco)~^, which is usually 
much smaller than the boundary length L. Hence, although this is strictly a q = instability, one expects in 
practice to observe a densely wiggling boundary on the scale of the entire domain. Another interesting observation is 
that, beyond the onset of rippling, q* does not always increase monotonically with compression. The nonmonotonic 
behavior becomes more pronounced the smaller the value of t/{KScq), as demonstrated in Fig. ^. For small values of 
this parameter, therefore, one expects the boundary to ripple densely beyond the onset of instability and then, upon 
further compression, return to a less rough shape. Recall that our Ansatz concerning the preferred perturbation is 
expected to give reliable results as long as <? ^ A^^. We have found that the ripphng modes obey q < Scq. On the 
other hand, a stable, sharp topography requires A^^ Scq/A (see Eq. (||)). Thus, our assumption is only marginally 
fulfilled and the results should be regarded merely as a qualitative guide. 
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VI. BEHAVIOR NEAR A CRITICAL POINT 



The topographic effects described in this article rely on a contrast between different domains. Hence, when a 
monolayer at coexistence reaches a critical point, these effects are expected to vanish along with the domain structure. 
Various parameters affecting the topography dramatically change when the critical point is approached: the density 
contrast becomes increasingly weak (leading to a smaller 6co), domain boundaries get diffuse (larger ^), and the bare 
line tension between domains, r, tends to zero. Thus, although the topography must clearly disappear at the critical 
point, its exact behavior has to be examined in detail. For example, it is unclear a priori whether, with respect to 
topography, the monolayer is driven towards the diffuse limit (larger ^Scq) or the sharp one (smaller ^Scq). 

As the temperature T approaches its critical value Tc, T = \T — Tc\/Tc ^ 0, we have fs^ 

Sca - f '3 0, (30) 
where, for a two-dimensional fluid, 1^ = 1 and f] — 1/8. Hence, the height difference, given by Eq. (|l|), decays as 

h = X^Sco r^f^^ f (31) 



Since S^Scq ~ j^-i^+l^ ^ y'-7/8 _^ j-j- jg ^^^c diffuse limit of Sec. IV that applies near the critical point. (The three- 
dimensional topography could affect the critical behavior of the two-dimensional ffuid in the monolayer as a "hidden" , 
annealed variable. Hence, the critical exponents should be modified according to the Fisher renormalization [ p3[ . In 
the case of a two-dimensional fluid (or Ising model), however, the Fisher renormalization leaves the exponents intact.) 

We now explore further details of the topographic critical behavior. For small inflection angles we expect in the 
diffuse hmit (cf. Eq. (||)) 

00 ~ h/^ - e+l^ ^ f (32) 



Indeed, the topography has been found in Sec. IV to depend on two dimensionless quantities, £,Scq and ^/{X^Scq) — ^/h, 
both of which diverge at the critical point — the former as T^^^^ — T^'^/^ and the latter as T~'^~^ = T-^/^, Studying 
Eqs. ( p^ ) and ( |20| ) in this asymptotic limit, one finds 6'o — /i/C and 9i ~ 0o/2, which verifies Eq. (32). Substituting 
these results in Eq. (|l8|) for the infiection energy, we get 
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Thus, the topographic contribution to the heat capacity of the monolayer vanishes as dg/dT ^ T^^^, whereas the heat 
capacity of the two-dimensional fluid diverges logarithmically ||3^ . This consistently demonstrates that the critical 
behavior of the monolayer remains unaffected by the topography. 

How does the approach to the critical point influence the instabilities studied in the previous sections? The folding 
instability in the diffuse limit requires, according to Eq. (p^), a surface tension lower than 



7e ~ Kdco/^ - f''+^ - f ^ 0. (34) 

In practice, therefore, as soon as the required lateral pressure exceeds that of the critical point the folding instability 
will become unattainable. 

The effect on the boundary-rippling instability studied in Sec. ^ is more delicate. Upon approaching the critical 
point the bare line tension of the boundary gets vanishingly small as [ ^4|j35[] 



T - 0, fi^l. (35) 

Thus, the resistance to rippling becomes increasingly weak. Yet, at the same time the driving force for rippling, i.e., 
the energy g gained due to the inflected topography, gets weaker as well. According to Eq. (^3|) the latter vanishes 
slightly faster, as T^^/^, Hence, it is the bare line tension that wins close to the critical point, and the boundary 
topography should flatten out at Tc as a smooth step without ripples. 
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VII. DISCUSSION 



We have demonstrated in this work that biphasic monolayers are generally nonflat, having inflected shapes in the 
vicinity of domain boundaries. This leads to an overall topography of mesas where domains of one phase are higher 
than those of the other. As the monolayer is progressively compressed the mesas grow more pronounced, subsequently 
developing overhangs, and finally becoming unstable. 

Substituting typical values for phosphohpid monolayers Q — 7 ~ 10-50 erg/cm^, K ~ 10-50 ksT, Cq ^ ~ 5-10 
nm — we get A ~ 1-10 nm, XScq ~ 0.1-1, and h ~ 0.1-10 nm. Hence, the mesas are steep but low. The numerical 
value of XScq implies that the predicted instability {XScq > 4) may be observed for attainable pressures. The energy 
per unit length gained by departing from the flat state to a sharp inflection is, according to Eq. (H), g ~ KScq ~ 1-10 
k^T /nm. Hence, for a typical domain size of L ~ 1-10 /im, the inflected conformation is "frozen", i.e., robust under 
thermal fluctuations. This justifies our mechanical, "zero-temperature" approach. 

As a more specific example, we may consider a monolayer whose behavior is governed by electrostatic interactions. 
The deviation from a fiat conformation is thereby driven by variations in the lateral charge density, cr(s). In the 
typical case of strong screening, c/k <C 1, where is the Debye screening length, one obtains ||] Kcq = 7rcr^/eK^ 
and K = Sira'^ /2eK^ , e being the dielectric constant of water. (Note the finite, positive Cq; charged monolayers 
spontaneously tend to curve into the aqueous phase.) Consequently, substituting typical values of a ~ 1 charge per 
0.3-1 nm^ and ~ 1-10 nm, we reach similar conclusions to those above. 

We have studied domain boundaries of finite thickness as well. The qualitative features of infiected conformation 
and instability do not disappear for any boundary thickness ^. On one hand, for a given compression a diffuse 
boundary leads to more moderate slopes compared to a sharp one. On the other hand, it shifts the folding instability 
to a higher pressure, thus strengthening the monolayer and allowing for higher mesas to be stabilized. Unfortunately, 
conventional means of increasing ^, e.g., heating towards a critical point, also reduce the domain contrast Sco, thus 
suppressing the topography. We have studied this delicate interplay close to a critical point in Sec. VI. Our results 
for diffuse boundaries show that the simple, infinitely sharp limit gives good results as long as ^ < (Scq)^^, which 
holds in most practical circumstances except near a critical point. 

One might worry about additional factors that would destroy the inferred topography. Such a factor is the cost 
in gravitational energy of displacing water from the flat interface. This energy per unit area is about 5pgh? ~ lO'' 
k^T/cTT? , where 5p is the difference in density of the two phases and g is here the gravitational acceleration. Thus, 
due to the small height of the mesas (1-10 nm), gravity is negligible over all relevant lateral length scales (up to 
meters). (Beyond the topographic instability, however, the monolayer may become much more folded, and gravity 
may have a significant stabilizing role.) Another factor to worry about is the van der Waals attraction between the 
inferred overhang and the underlying surface, which might make the overhang collapse. The attraction energy per 
unit area is roughly Hjh? , where H is the Hamaker constant divided by 127r (typically a few k-gT) and h ~ X^Scq is 
the overhang height g^l- The lateral extent of the overhang is A, and the resulting energy per unit length, H/ {X^6c^), 
is to be compared with the infiection energy, K6cq. The ratio is [H / K){X5c[))~^ <C 1, since K oi & lipid monolayer is 
a few tens k^T and X5co ~ 3-4 to get an overhang (cf. Eq. (||)). Hence, the van der Waals attraction is too weak to 
significantly affect the overhang. 

The topography of mesas and overhangs is thus a robust result which should be observable in practice. Such an 
observation is difficult, however, because of the small height differences and fiuidity of the interface. Very recently a 
new experimental technique has been presented, utilizing non-specular scattering of intense light to visualize small 
topographic features in phopholipid monolayers js^. Although the study was focused on features of a pure liquid- 
condensed phase, height differences were reported at boundaries of liquid-expanded domains coexisting with a gas 
phase, as well as liquid-condensed domains in a liquid-expanded phase. (Interestingly, a stronger signal was obtained 
in the former case, perhaps due to a larger contrast in spontaneous curvature.) It is still unclear whether these 
experimental findings are related to the topography discussed here or to other, more molecular effects. 

Recent experiments on mixed phospholipid monolayers have revealed a new type of folding instability pl| , p3[ . When 
the monolayer is compressed and enters a coexistence region, there is a critical pressure at which micron-scale folds 
appear. The folding is significantly more reversible than other collapse mechanisms and is therefore thought to be of 
key importance to the function of lungs. Figure |^ shows a pressure-area isotherm as measured for a mixed phospholipid 
monolayer of dipalmitoylphosphatidylcholine (DPPC) and palmitoyloleoylphosphatidylglycerol (POPG). The folding 
is manifested by a plateau in the isotherm. (The same phenomenon was observed in DPPG monolayers at a much 
higher surface tension ||2l[].) Figure ^ presents a sequence of fiuorescence microscopy images of the monolayer just 
before and just after the instability. 

We believe that this folding phenomenon is initiated by the topographic instability of boundary regions as obtained 
from our model. (Further evolution and propagation of the fold are determined by other factors not taken into 
account in the current work, such as the viscoelasticity of the monolayer [§3|,58l.) If the hypotheses regarding the 
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biological significance of the folding and its relation to topography are correct, it may represent an interesting solution 
of Nature to a delicate mechanical problem. Using a mixed surfactant monolayer to cover the lung leads to domain 
formation upon compression, which in turn allows the topographic instability and folding. Additional constituents 
{e.g., proteins) may ensure that the folding is not preceded by other, irreversible collapse mechanisms This 
design provides the monolayer with a unique way to yield gracefully to compression and reduce its projected area, 
while avoiding irreversibility and loss of surfactant. 

Folding of the mesa structure is in many cases preempted by other instabilities. One type of collapse is 
delamination — breakage of the monolayer into multiple layers [pyiO|. It occurs when the surfactant sheet yields 
to a combination of bending and lateral compressive stresses. Since the mesas help relieve part of the inherent bend- 
ing stresses exerted in a flat monolayer, one expects the breakage to occur (somewhat counter-intuitively) away from 
the boundary, inside the more frustrated domain (i.e., the one having higher spontaneous curvature). Another mode 
of monolayer collapse is budding and ejection of vesicles into the aqueous phase |Tl| . Recent experiments on mixed 
phospholipid monolayers have shown that vesiculation is promoted by increased temperature and may coexist with 
folded structures The effect of topography on delamination and budding, as well as the interplay between the 
various collapse modes, are yet to be studied in detail. 

Another general conclusion arising from this work relates to fluid surfaces of vanishing tension. Such surfactant films 
are encountered, e.g., in emulsions, L3 ("sponge") phases and large, unsupported bilayer vesicles The topographic 
instability found for a finite tension implies that these tensionless surfaces cannot sustain a stable domain structure. 
To the best of our knowledge a static domain structure has never been observed in those systems. We attribute it to 
the inevitable shape instability that would occur near domain boundaries if such a structure existed. 

We have studied a possible rippling of the domain boundary upon compression. This phenomenon arises from a 
competition between the topographic features accompanying the domain boundary and its bare line tension. The 
threshold value of line tension required to get rippling is Tc ~ ~ 1-10 fceT'/nm, which is of the same order as 
line tension values measured in experiments |40-44j. The rippling is therefore a realistic, observable feature. There 
is already a well-established mechanism for shape transformations of monolayer domains, driven by a competition 
between line tension and long-range electrostatic interactions [^5|" ^7|j40| . We offer the topographic rippling as an 
additional mechanism which should be observed in practice. There are three major features distinguishing the two 
phenomena, (i) In the electrostatic mechanism an infinite straight boundary is never stable. Consequently, the stable 
domain size and wavelength of boundary instability have the same scale, L ^ ^ le^/^P , where / is a molecular 
size and 5p the difference in dipole densities of the two phases [^. Hence, shape deformations occur on the scale of 
the entire domain [HO 
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4q ], leading to a sequence of well-resolved, "quantized" transitions. By contrast, the length 
scale of the topographic rippling, (9*)"^ ~ (5ca)~^ ^ 0.1 /im, is much smaller than, and unrelated to, L. Thus, we 
expect this instability to appear as a small-scale roughening of the domain boundary, (ii) The topographic rippling, 
being an elastic mechanism, should not be very sensitive to changes in electrostatic parameters such as ionic strength 
and molecular charge, (iii) As demonstrated in Sec. 0, the rippling wavelength may exhibit in certain circumstances 
a peculiar nonmonotonic behavior as a function of pressure. In addition, we have shown in Sec. VI that topographic 
rippling is inhibited near a critical point. Hence, boundary rippling could be smoothed out by heating the monolayer 
towards its critical temperature. This is in contrast with common surface instabilities that are usually promoted by 
increasing temperature. 

In a recent experiment on a pure DPPC monolayer at liquid expanded-liquid condensed coexistence, a domain- 
boundary instability of sub-micron scale has been observed [p2| . This is demonstrated in Fig. ^ At a critical pressure 
slight roughening appears simultaneously in all domain boundaries. Upon little further compression the roughening 
becomes denser and the boundaries look fuzzy due to optical limitations. A detailed presentation of this effect will 
be given in a forthcoming publication . The small length scale of this instability (compared to the entire domain 
size) is in accord with the topography-induced mechanism discussed above. Yet, further study is required in order to 
clarify the relation between the two effects. 

The phenomena described in this work — mesa formation, mesa instability, boundary rippling — arise from rather 
basic considerations. Nevertheless, there is still a gap between theory and experimental observations. The relation 
between the topographic instability as obtained from the elastic model and the observed folding in biphasic lipid 
monolayers is still to be established. In particular, the current theory does not account for the fully developed folded 
structure and its stability, as observed in experiments. Topography-induced boundary rippling and its distinction 
from the known electrostatic mechanism is yet another intriguing feature to be experimentally investigated. We hope 
to close this gap in future publications. 

The mesa topography is a novel interfacial feature predicted by our work. If mesas exist, which is yet to be 
decisively proven by experiment, they should have important implications on various interfacial aspects, such as surface 
interactions with dissolved molecules, behavior in confined geometries, and possible applications for controllable 
nanostructures. We would like to draw special attention to the unique overhang topography predicted by the model. 
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Under the right compression all the domains in a biphasic monolayer should develop regular lips at their edges. Such 
controllable nanoscale grooves might be technologically useful, e.g., for capturing and encapsulating (bio)polymers. 
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APPENDIX: DIFFERENTIAL GEOMETRY OF A RIPPLED BOUNDARY 

In the rippled state we represent the domain boundary by the curve R(i/) = {x = a sin qy, y, z = 0) and parameterize 
the monolayer surface as 

r{s,t) = {x{s,t),y{s,t),z{s,t)) 

x(s,t) ^ / cos 9 {s' )ds' + a sin qt, y(s,t) = t, z{s,t) ^ / sin 6'(s')ds'. (Al) 



We now need to represent the various properties of the surface using the " mate rial coordinates" {s,t) [Q. 



The determinant of the metric tensor associated with the surface of Eq. (Al) is 



F = (dsT X dtv)^ = 1 + {qa cos qt sine f. (A2) 

The Jacobian of the transformation {x,y) is,t) is J ~ dgxdty — dixd^y = cos9. Using these expressions we find 
the area element, 

d^(s, t) = F^/^dsdt = (1 + q^a^ cos^ qt sin^ Of/'^Asdt, (A3) 
and its projection onto the xy plane, 

dAp(s, t) = Jdsdt = cosOdsdt. (A4) 
An element of the boundary curve is given by dR = {qa cos qt, 1, 0)dt, and the resulting arc- length element is 

di?(s,t) = [l + {qacosqtf]^/'^dt. (A5) 
What is left to calculate is the local surface curvature, c{s,t). The local normal to the surface is given by 

n = F"^/^(a,r X dtr) = r^^/^(- sine,qacosqtsin0,cose). 
The mean curvature can then be calculated either from the trace of the curvature tensor, 



c{s,t) = --tT{dni/drj), 



or by momentarily resorting to the Monge representation, 

c^V ■[T-'^^Vz{x,y)], 
where the V operator is defined in the xy plane. The result is 

c(s, t) = r~3/2 [(1 + g^a^ cos^ qt)9 + q^a sin qt sin 9] . (A6) 
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FIG. 1. A) Schematic sketch of the monolayer in the boundary region. A boundary lying parallel to the y axis separates two 
large domains denoted by 1 and 2. B) Cross-section parallel to the xz plane. The monolayer conformation is parameterized by 
the angle 9{s) it makes with the xy reference plane at curvilinear distance s from the boundary. 
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FIG. 2. A) Slope angle profiles near a sharp domain boundary as compression is increased. The curves are obtained from 
Eqs. (^) and (^ using the values (from bottom to top) \5cq = 1, 3, 4. B) The corresponding spatial conformations. All lengths 
are given in units of {5cq)~^ . 
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FIG. 3. Effect of diffuse boundary on mesa instability. A) Rescaled critical tension, (^/i4'5co)7c, required for tlie transition. 
Tlie solid line is the result for a boundary of finite tliickness (i.e^tlie scaling function F{x) of Eq. (p3|)). Tlie dasfied line is 
the corresponding result for an infinitely sharp boundary (Eq. (P_2|)). B) Critical inflection angle. For very sharp boundaries 
the angle at instability is n, whereas for very diffuse ones it is reduced to n/2. 
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FIG. 4. Surface tension-line tension diagram of topographies for a sharp domain boundary. If r > Tc = (tt — 2)KSco the 
mesa topography remains straight upon decreasing surface tension until it becomes unstable at 7 = 7c = j^K{5co)'^ . If r < Tc 
the mesa wall ripples below a surface tension 7 = 7r > 7c- 
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FIG. 7. Prossurc-arca isotherm for a mixed monolayer of DPPC and POPG, as measured during a compression/expansion 
cycle in a Langmuir trough. The mole ratio is DPPC:POPG=7:3 and the temperature 25°C. The folding instability is indicated 
by an arrow. 




FIG. 8. Fluorescence microscopy images of the folding instability. A) Section of the monolayer just before folding {t = 0), 
exhibiting the biphasic domain structure. Dark regions are DPPC-rich; bright ones are POPG-rich. B) The same section 
at t =1/30 s. A micron-scale fold appears in between domain walls (indicated by arrow). The image is blurred because of 
monolayer movement during folding. The inset shows a contrast-enhanced image of the fold, magnified by 50 percent. C) The 
fold at t =2/30 s, having propagated to nearby domains. D) The fold at t =4/30 s, after the fast monolayer movement has 
ceased. 
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FIG. 9. Roughening of domain boundaries upon compression, as observed by fluorescence microscopy. The monolayer consists 
of DPPC phosphohpids and hcs over au aqueous solution of 0.2M NaCl at temperature 24.5° C. Dark regions correspond to a 
liquid-condensed phase and bright ones to a liquid-expanded phase. The lateral pressures and molecular areas (in mN/m and 
A^, respectively) are: A) 16.9, 54.8; B) 19.4, 53.0; C) 21.6, 51.9; D) 23.4, 51.2. Images were contrast-enhanced. 
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